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ABSTRACT 

Context. Asteroid (21) Lutetia is the target of the ESA Rosetta mission flyby in 2010 July. 

Aims. We seek the best size estimates of the asteroid, the direction of its spin axis, and its bulk density, assuming its 

shape is well described by a smooth featureless triaxial ellipsoid, and to evaluate the deviations from this assumption. 

Methods. We derive these quantities from the outlines of the asteroid in 307 images of its resolved apparent disk obtained 

with adaptive optics (AO) at Keck II and VLT, and combine these with recent mass determinations to estimate a bulk 

density. 

Results. Our best triaxial ellipsoid diameters for Lutetia, based on our AO images alone, are a x b x c — 132 x 101 x 93 

km, with uncertainties of 4 x 3 x 13 km including estimated systematics, with a rotational pole within 5° of ECJ2000 

[A P] = [45° - 7°], or EQJ2000 [RA DEC] = [44° + 9°]. The AO model fit itself has internal precisions of 1 x 1 x 8 

km, but it is evident, both from this model derived from limited viewing aspects and the radius vector model given 

in a companion paper, that Lutetia has significant departures from an idealized ellipsoid. In particular, the long axis 

may be overestimated from the AO images alone by about 10 km. Therefore, we combine the best aspects of the radius 

vector and ellipsoid model into a hybrid ellipsoid model, as our final result, of 124 ± 5 x 101 ± 4 x 93 ± 13 km that can 

be used to estimate volumes, sizes, and projected areas. The adopted pole position is within 5° of [A /3] — [52° — 6°] or 

[RADEC] = [52° -)-12°]. 

Conclusions. Using two separately determined masses and the volume of our hybrid model, we estimate a density 

of 3.5 ± 1.1 or 4.3 ± 0.8 g cm"''. From the density evidence alone, we argue that this favors an enstatite-chondrite 

composition, although other compositions are formally allowed at the extremes (low-porosity CV/CO carbonaceous 

chondrite or high-porosity metallic). We discuss this in the context of other evidence. 
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■^ 1. Introduction 



The second target of the ESA Rosetta mission, asteroid 
(21) Lutetia, had a favorable opposition in 2008-09, reach- 
ing a minimum solar phase angle of cj = 0.45° on 2008 
November 30, and a minimum distance from the Earth of 
1.43 AU a week earlier. Based on previously determined 
sizes of Lutetia, with diameters ranging from 96 km from 
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IRAS [Tedesco et al. 2002, 2004] to 116 km from radar 
Magri et al. [1999, 2007], Lutetia should have presented an 
apparent diameter of 0.10", slightly more than twice the 
diffraction limit of the Keck Observatory 10 m telescope at 
infrared wavelengths (1-2 /im). Continuing our campaign 
[Conrad et al. 2007; Drummond et al. 2009; Carry et al. 
2010a] to study asteroids resolved with the world's large 
telescopes equipped with adaptive optics (AO), we have 
acquired more than 300 images of Lutetia, most from the 
2008-09 season. An exceptionally good set of 81 images was 
obtained on 2008 December 2 with the Keck II telescope, 
which, despite the high sub-Earth latitude, yields a full tri- 
axial ellipsoid solution from the changing apparent ellipses 
projected on the plane of the sky by the asteroid. Analyzing 
all available images (2000, 2007, and 2008-09 seasons) yields 
a result consistent with the December 2 set. 
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Fig. 1. Sample AO images of Lutetia at nine epoclis on 
2008 December 2. The scale is in pixels, where 1 pixel cor- 
responds to 0.01". North is up and east is to the left. Images 
are displayed on a linear scale, after bi-linear interpolation. 
From left to right, top to bottom, the images were obtained 
between 6:52 and 9:01 UT after subtracting for light-time 
travel. They correspond to the 9 points in Fig 3. 



2. Observations 

Table 1 gives the observing circumstances for all seven 
observation dates, where, in addition to the date. Right 
Ascension, Declination, and Ecliptic longitudes and lati- 
tudes of Lutetia, we give its Earth and Sun distance, the 
solar phase angle (w), the position angle of the Sun mea- 
sured east from north while looking at the asteroid (NTS), 
the filter and the number of images (to) used to form mean 
measurements at n epochs, the scale in km/pix at the dis- 
tance of the asteroid on the date, and the instrument and 
telescope configurations used for the observations as listed 
in Table 2. Multiplying the first scale in Table 1 by the 
scale in Table 2 will give a km/" scale at the distance of 
the asteroid, and multiplying this by the appropriate res- 
olution element (Q) from Table 2, which is the diffraction 
limit = X/D, where A is the wavelength and D is the 
telescope diameter, gives the last column of Table 1, the 
km per resolution element scale. 

The best data set, obtained at K' on December 2, com- 
prises images at 9 epochs, 9 images per epoch, where each 
image is a 0.4 second exposure. Figure 1 shows a single 
image from each of the 9 epochs, and clearly reveals the 
asteroid rotating over a quarter of its 8.2 h period. 

Figure 2 illustrates the range of solar phase angles for 
our observations, showing that on more than half of the 
dates the phase angle was greater than 17°, which, for ir- 
regular bodies, can lead to violations of the assumptions 
we adopt in Section 3. In this particular case, however, we 
have found the data are still quite useful and contribute 
substantially to our final results. 
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Fig. 2. Range of solar phase angles (w) for (21) Lutetia 
observations. All (circles) but two dates (squares) are from 
the 2008-09 opposition. We list the date and observatory for 
observations obtained outside the 2008-09 campaign. The 
number of epochs for each date are listed in parentheses. 



3. Analysis 

Assuming that an asteroid can be modeled as a uniformly 
illuminated triaxial ellipsoid rotating about its short axis, 
it is possible to estimate its diameters {a > b > c) and find 
the direction of its spin axis from the observation of a series 
of ellipses projected as it rotates. We can thus treat any as- 
teroid in the same manner as we treat asteroids that are 
well described by the triaxial assumption, e.g., those that 
we have defined as Standard Triaxial Ellipsoid Asteroids, or 
STEAs [Drummond ct al. 1985, 1998, 2009; Conrad et al. 
2007; Drummond & Christou 2008]. The key to turning the 
projected ellipses into a triaxial ellipsoid is determining the 
ellipse parameters from AO images. We use the method of 
Parametric Blind Deconvolution [PBD, Drummond et al. 
1998; Drummond 2000] to find the long (a) and short (/3) 
projected (plane of sky) ellipse axes dimensions and the 
orientation or position angle (PA) of the long axis. PBD 
allows us to find not only the asteroid ellipse parameters 
but parameters for the Point Spread Function (PSF) as 
well. Having shown that a good model for the AO PSF is 
a Lorentzian [Drummond et al. 1998; Drummond 2000] we 
simply treat a disk-resolved image of an asteroid as the con- 
volution of a flat-topped ellipse and a Lorentzian, making 
a simultaneous fit for each in the Fourier plane where the 
convolution becomes a simple product. 

All 307 images were fit for the projected asteroid ellipse 
parameters and individual Lorentzian PSFs. The mean and 
standard deviations of the parameters were formed around 
epochs consisting of a series of 3 to 15 images obtained in 
one sitting at the telescope. We then solve the triaxial ellip- 
soid from a least square inversion of the ellipse parameters 
[Drummond et al. 1985; Drummond 2000]. 

The projected ellipse parameters can also be extracted 
from images deconvolved with an alternate algorithm such 
as Mistral [Conan et al. 2000; Mugnier et al. 2004]. Such 
contours provide a direct measurement of details of the 
projected shape of the asteroid, allowing the construction 
of the radius vector model that we present in Carry et al. 
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Table 1. Observation Log for (21) Lutetia 



Date 


EQJ2000 


ECJ2000 


V 


Earth 


Sun 


OJ 


NTS 


Filter 


Scalel 


Config 


Scale2 


(UT) 


(RA° Dec°) 


(Lon° Lat°) 




(AU) 


(AU) 


n 


n 


(m/n) 


(km/pix) 


Table 2 


(km/rex) 


Aug 15, 2000 


12.3 - 0.3 


11.2 -5.1 


10.5 


1.239 


2.057 


21.6 


71.6 


K'(19/l) 


15.10 


A 


39.8 


Jun 6, 2007 


246.0 -20.7 


247.8 +0.9 


10.1 


1.295 


2.305 


3.3 


286.4 


Ks(35/7) 


12.47 


B 


51.7 


Oct 22, 2008 


75.1 +20.9 


76.2 -1.8 


11.2 


1.555 


2.352 


17.9 


85.9 


J (15/2) 


11.21 


C 


29.3 


Oct 22, 2008 


75.1 +20.9 


76.2 -1.8 


11.2 


1.555 


2.352 


17.9 


85.9 


H (15/2) 


11.21 


C 


38.3 


Oct 22, 2008 


75.1 +20.9 


76.2 -1.8 


11.2 


1.555 


2.352 


17.9 


85.9 


K'(33/3) 


11.21 


C 


49.6 


Nov 21, 2008 


69.4 +20.8 


70.9 -1.3 


10.5 


1.430 


2.406 


4.7 


88.5 


K'( 4/1) 


10.31 


C 


45.6 


Dec 2, 2008 


66.4 +20.6 


68.1 -1.1 


10.2 


1.441 


2.426 


1.1 


236.7 


K'(81/9) 


10.39 


C 


46.0 


Jan 23, 2009 


59.6 +20.6 


61.9 +0.1 


11.8 


1.895 


2.518 


20.1 


258.5 


K'(30/2) 


13.66 


C 


60.5 


Feb 2, 2009 


60.6 +20.9 


62.8 +0.2 


12.0 


2.033 


2.534 


21.5 


258.9 


H (30/2) 


14.66 


C 


81.1 


Feb 2, 2009 


60.6 +20.9 


62.8 +0.2 


12.0 


2.033 


2.534 


21.5 


258.9 


K'(45/3) 


14.66 


C 


64.9 



Table 2. Configuration for Table 1 



Configuration 


Instrument 


Telescope 


Aperture 


Scale 


Filter 


Wavelength 


Resolution 








(m) 


(Pix/") 




(^m) 


(") 


A 


NIRSPEC 


Keck II 


10 


59.52 


K' 


2.12 


0.044 


B 


NACO 


ESO VLT UT4 


8.2 


75.36 


Ks 


2.18 


0.055 


C 


NIRC2 


Keck II 


10 


100.6 


J 
H 
K' 


1.25 
1.63 
2.12 


0.026 
0.034 
0.044 



[2010b], providing a more refined description of the shape 
of Lutetia. However, only the PBD parameters were used 
in deriving ellipsoid solutions here. 

While observations from one night can pro- 
duce triaxial ellipsoid results [Conrad et al. 2007; 
Drummond & Christou 2008; Drummond et al. 2009] 
it is possible to combine observations from different nights 
over multiple oppositions to make a global fit if a sidereal 
period is known with enough accuracy (Drummond et 
al. in preparation). This can resolve the natural two- fold 
ambiguity in the location of a rotational pole from a single 
night of data (see Section 4.2), and in some cases reduce 
the uncertainty in the triaxial ellipsoid dimensions if the 
asteroid is observed over a span of sub-Earth latitudes. 
For instance, when observations are restricted to high 
sub-Earth latitudes, even a few images at an equatorial 
aspect will supply a much better view of the c axis than a 
long series at near polar aspects. In other words, different 
viewing geometries generally lead to a better solution. 

Unfortunately, during the last two oppositions, in 2007 
and 2008-09, Lutetia's positions were ^ 180° apart on the 
celestial sphere. In fact, the position of Lutetia for the 
Very Large Telescope (VLT) observations on 2007 Jun 6 
was exactly 180° from its position on 2008 December 2 
(see Table 1), which meant that regardless of the location 
of the rotational pole, the two sets of observations were 
obtained at the same sub-Earth latitude, but of opposite 
signs. Thus, with the now-known pole, our observations 
of Lutetia over the 2008-09 opposition were at the same 
(but southerly) deep sub-Earth latitudes as the VLT ob- 
servations in the previous (but northerly) deep sub-Earth 
latitude opposition in 2007, affording the same view of 
the strongly fore-shortened c axes, but from the opposite 
hemisphere. The single set of images in 2000 was also 
obtained at deep southerly latitudes. Thus, we have no 
equatorial view of Lutetia and, therefore, its shortest (c) 
dimension remains less well determined than the other 
two. 



Table 3. (21) Lutetia Ellipsoid Fit Solutions 





Triaxial (Dec 2008) 


Biaxial (All) 


a (km) 


132+1 


132+1 


b (km) 


101+1 


101+1 


c (km) 


93 +8 


c=b 


en 


-66±3 


-59 + 3 


PANodeC) 


155+3 


178+6 


^0 (Max) (UT) 


5.66+0.07 


5.10+0.12 


Pole 






[RA°; Dec°] 


[44;+9] 


[34; +16] 


a radius (°) 


2.5 


3.1 


[A°;/3°] 


[45;-7] 


[37; +3] 



4. Results 

4.1. Ellipsoid Fits 

We made two separate fits to our data, one using only our 
best data set from December 2, and the other combining 
all of our data taken over the 8.5 year span using a sidereal 
period of 8.168 27 h from Carry et al. [2010b] to hnk them 
together. The results of these two fits give highly consistent 
values for the equatorial dimensions, but not for c. When 
fitting all of the data, a higher c value is preferred, but we 
restrict c to the usual physical constraint of & > c. Because 
the entire data set, taken as an ensemble, has higher noise, 
we consider the 6 = c as a limiting case (a biaxial ellipsoid), 
and adopt the value of c derived from the Dec 2 data alone 
as our preferred value. Both solutions are listed in Table 3, 
where 9 is the sub-Earth latitude, PAjsiode is the position 
angle of the line of nodes measured east from north, and 
ipo is rotational phase zero, the time of maximum projected 
area when the a axis lies unprojected in the plane of the sky. 
The uncertainties shown in Table 3 are the internal preci- 
sions of the model fit, and do not include systematics. We 
have assigned overall uncertainties to our ellipsoid model, 
including systematics, of 4 x 3 x 13 km in the dimensions, 
and 5° in the pole position. 
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Table 4. RMS of Projected Ellipses from Models 



2008 Dec 2 UT sub-Earth Lat— 66 Solar phase-1 .1 



Triaxial (2008 Dec only) 
Triaxial (All) 
Biaxial (All) 



Of 

(km) 



0.4 
4.1 
3.4 



/3 
(km) 



0.9 
4.7 
4.4 



PA 

n 



1.1 

6.9 
6.6 



Fie 



Figure 3 shows the triaxial ellipsoid fit to the December 
data. Figure 4 shows the residuals from the prediction us- 
ing the December triaxial ellipsoid model for all of the data, 
and Figure -5 shows the overall residuals to the biaxial fit. 
The difference between Figs 4 and 5 is subtle, showing that 
Lutetia is close to a prolate ellipsoid. Trends in some of the 
residuals, such as the curling set of position angles at the 
right of the top plot in Fig 4, and at the left of the simi- 
lar plot in Fig 5, indicate departures from our assumptions 
of a smooth featureless ellipsoid rotating about its short 
axis, and motivate the more detailed shape model that we 
present in Carry et al. [2010b]. The rms weighted (by the 
observational uncertainty of each measurement) residuals 
for a and (3, the projected ellipse major and minor axes 
lengths, and the weighted residuals for the position angle 
(PA) of the long axis, are given in Table 4, and are to be 
associated with Figs 3-5. These can be interpreted as un- 
certainties (but without possible systematics) for any pre- 
dicted future projected ellipse parameters. 

The axial ratios derived from our model are a/b =1.32 
and b/c — 1.09. From a compilation' of axial ratios and ro- 
tational poles, mostly from lightcurves, the average axial ra- 
tios are a/b = 1.27±0.06 and b/c = 1.45±0.55, both within 
one sigma of our directly determined values. More recent 
work suggests a b/c ratio of less than 1.1 [Belskaya ct al. 
2010], also consistent with our value. Our fit of the ra- 
dius vector model derived from lightcurves by Torppa et al. 
[2003, see section ^) below] yields an a/b ratio of 1.17 and 
b/c = 1.18, and the latest radius vector model derived by 
Carry et al. [2010b] from a combination of present AO im- 
ages and lightcurves has ratios of 1.23 and 1.26. Our hybrid 
model, given as the final result of our paper here (see next 
section), has a/b = 1.23 and b/c = 1.09. 

4.2. Rotational Pole 

There is a natural two-fold ambiguity in the location of the 
rotational pole with our method that is symmetric with re- 
spect to the position of the asteroid if observed on one night. 
Thus, there are two possible poles for the December triaxial 
ellipsoid solution. However, if the asteroid can be observed 
at significantly different positions, then the rotational pole 
can be disambiguated. The single 2000 observation helps 
break the ambiguity since the residuals are some 18% higher 
for the rejected pole than for the accepted region when con- 
sidering all of the data. Otherwise, the 2008-09 observations 
and the 2007 data from the same positions, would not have 
provided enough diversity to break the ambiguity. 

The poles from various lightcurve techniques can have 
two- or four-fold ambiguities [see Magnusson et al. 1989, 
for a good summary], which can be broken when paired 
with our results. Figure 6 shows the positions of about 



^ http://vesta.astro.amu.edu.pl/Scicncc/Astcroids/ is a web 
site gathering sidereal periods, rotational poles, and axial ratios 
maintained by A. Kryszczyriska. See Kryszczynska et al. [2007]. 




Rotational Phase 



Fig. 3. Triaxial ellipsoid fit to measured ellipse parameters 
of (21) Lutetia on 2008 Dec 2. In the upper subplot, each 
image's long (a) and short (/?) axis dimensions are plotted 
as dots with 1 a uncertainties. The lines are the prediction 
for the projected ellipses from the triaxial ellipsoid parame- 
ters in Table 3, derived from the fit to the data. Because the 
solar phase angle (w) is only 1.1° (Fig. 2 and Table 1), the 
ellipse parameters for the terminator ellipse is coincident 
with the projected ellipse lines. The lower subplot shows 
the same for the position angle of the long axis, where the 
horizontal line is the line of nodes, the intersection of the 
asteroid's equator and the plane of the sky. This figure is 
corrected for light time travel, i.e., the plot is in the body- 
centered time frame. 
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Fig. 4. Residuals between the triaxial ellipsoid model (ob- 
tained from December 2008 observations only) and all ob- 
servations as function of the rotation phase and latitude. 
The major and minor axes residuals are shown as squares 
and circles, respectively, both in km, and the position angle 
residuals, in degrees, are shown as diamonds. From Table 4, 
the rms of the weighted residuals are 4.1 and 4.7 km for the 
major and minor axes, respectively, and 6.9° for the posi- 
tion angle. The data in the northern hemisphere is from the 
VLT in 2007. 
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Fig. 5. Same as previous figure, except that the residuals 
are between the biaxial ellipsoid model and the data. From 
Table 4, the rms of the weighted residuals are 3.4 and 4.4 
km for the major and minor axes, respectively, and 6.6° for 
the position angle. 



half of the poles (see footnote 1) found from lightcurve 
methods (the other half lie on the opposite hemisphere), 
as well as ours. We assert that the correct region for the 
pole location is where our poles near RA = 45° coin- 
cide with the span of lightcurve poles in this hemisphere. 
Furthermore, the lightcurve inversion (LCI; see section 5) 
pole of Torppa et al. [2003] " lies at coordinates RA= 52° 
and Dec= +13°, less than 9° from our triaxial ellipsoid pole 
in Table 3. 



5. Comparison with Lightcurve Inversion IVIodel 

Figure 7 shows our PBD images of Lutetia from 2008 
December 2. Each image is the mean of 9 shifted and added 
images at each epoch, and then linearly deconvolved of the 
Lorentzian PSF found in its fit. Notice the tapered end. 

For comparison, the lightcurve inversion model (see 
footnote 2) based on the work of Torppa et al. [2003] 
is shown in Fig 8 for the same times. The model ap- 
pears to match the overall shape and orientation in the 
images, verifying the pole and sidereal period, but it 
does appear fatter and less tapered than the images 
of the asteroid. In a following article we combine our 
AO images with lightcurve data using a method known 
as KOALA (Knitted Occultation, Adaptive-optics, and 
Lightcurve Analysis [Carry et al. 2010b]) to produce an im- 
provement over the previous LCI model. Not only does it 
yield better matches to the AO images, but it provides 
an absolute kilometer scale, and it can reproduce Lutetia's 
lightcurve history. 

Although a triaxial ellipsoid fit of the new KOALA 
model yields diameters of 124 x 101 x 80, the model is 
very non-ellipsoid in appearance, and while the b diame- 
ter is in agreement between the AO-only and the KOALA 
model, both the KOALA a and c dimensions are ~ 10 km 
smaller than from our triaxial ellipsoid results here. The 




Fig. 6. Pole locations for Lutetia on the celestial globe. The 
positions of Lutetia on the sky for the seven nights of our 
observations are indicated by x 's. The opposite of the 2007 
position is the circle on top of the x for December 2008. The 
poles found from lightcurve work are marked with asterisks, 
with the LCI pole [Torppa et al. 2003] shown as a diamond 
(our final pole, the KOALA pole [Carry et al. 2010b], is less 
than a degree from the LCI pole). The four wedge shaped 
areas are the uncertainty regions around our poles, with T 
and B indicating our two each possible triaxial and biaxial 
ellipsoid solution poles. Our rejected poles are marked as 
Bx and Tx at the far right, while the lightcurve rejected 
pole region is on the other side of the globe. 
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Fig. 7. Same as Fig 1, except PBD (deconvolved) images 
of (21) Lutetia from 2008 December 2. 



http://astro.troja.mflF.cuni.cz/projects/asteroids3D/web.php 
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{abcY^^ = 105 ± 5 km. The rotational pole for these hy- 
brids should be the KOALA pole at [RA Dec] = [52° + 12°], 
since it is primarily based on lightcurves obtained over 47 
years as opposed to the triaxial ellipsoid pole in Table 3 
from nine epochs on 2008 December 2, although they are 
less than 8° apart. The uncertainty for this pole is about 
5° in each coordinate. 





o c o 



o 

Fig. 8. Convex shape model of (21) Lutetia from 
Torppa et al. [2003], plotted at same times as Fig 7. 



AO-only triaxial ellipsoid solution comes from only a quar- 
ter of a rotation, when the minimum area is projected (over 
what would be a lightcurve minimum). During the time of 
the 2008 Dec 2 AO observations, the b axis was seen un- 
projected in the plane of the sky, but both the a and c 
axes were not. It is the extrapolation, as it were, in ro- 
tation to the maximum projected area, when the a axis 
could be seen unprojected, that leads to an a dimension 
larger than found from KOALA. The KOALA technique, 
by combining lightcurves that cover all rotational phases 
and sub-Earth latitudes with the AO images (at restricted 
rotational phases and latitudes), finds that there is a large 
depression on the side of Lutetia away from the 2008 Dec 2 
observations that was not completely sampled by our imag- 
ing, resulting in the smaller KOALA model a axis dimen- 
sion. This depression of '^ 10 km explains the difference 
between the two a dimensions. 

On the other hand, since KOALA only uses amplitudes 
from lightcurves, and since amplitudes are a strong function 
of a/b but a weak function of b/c, the KOALA model c di- 
mension is only weakly determined when the AO images at 
high sub-Earth latitudes only reveal a strongly foreshort- 
ened c axis. (See [Carry et al. 2010b] for a discussion on 
the limits of the KOALA inversion in the particular case 
of Lutetia.) Therefore, Lutetia's c dimension is best deter- 
mined from the 2008 December 2 AO data set. 

To make the best possible model for use in evaluating 
sizes, cross-sectional areas, volumes, and densities, we com- 
bine what we consider the best aspects of both models into 
a hybrid triaxial ellipsoid/KOALA model that has dimen- 
sions of 124 X 101 X 93 km, taking the a diameter from 
KOALA and the c diameter from the AO triaxial ellipsoid 
fit. The original KOALA model radius vector Z components 
are merely expanded by 93/80. We estimate the uncertain- 
ties on these dimensions, including possible systematics, to 
be 5 X 4 X 13 km. Our best final average diameter is, then. 



6. Taxonomy and Density 

Lutetia was well observed in the 1970s, yielding visible and 
near-infrared reflectance spectra [McCord & Chapman 
1975], radiometric albedos and diameter estimates 
[Morrison 1977], and polarimetric albedos and diame- 
ter estimates [Zellner & Gradie 1976], which have been 
confirmed by similar observations reported during the 
last decade [see review by Belskaya et al. 2010]. Based 
on these data, Chapman et al. [1975] placed only three 
asteroids, (16) Psyche, (21) Lutetia, and (22) Kalliope, 
into a distinct taxonomic group to which Zellner & Gradie 
[1976] assigned the letter "M". The M type was defined 
in terms of spectral and albedo properties by Bowell et al. 
[1978], who assigned a diameter of 112 km to Lutetia 
(estimates by Morrison [1977] and Zellner & Gradie [1976] 
had been diameters of 108-109 and 110 km, respectively). 

It was later found by radar that some, but not all, 
M-types were metallic. Rivkin et al. [1995] recognized 
that there were two sub-types of M-type asteroids. The 
standard M types showed high radar refiectivity and 
relatively neutral colors, both apparently due to metal. 
The other type (also showing similar colors, but now 
thought to be due to metal flakes embedded in a colorless 
stony matrix) had a 3 micron band, ascribed to hydrated 
minerals and which was deemed to be unlikely on a chiefly 
metallic body. Rivkin et al. [1995] called this new "wet" 
subclass M(W) and assigned Lutetia to this subclass 
[Rivkin et al. 2000]. Chapman & Salisbury [1973] first 
suggested that what we now term an M-type spectrum 
might be associated with enstatite chondrites (ECs) and 
Rivkin et al. [2000] suggested a hydrated EC as a plausible 
composition for Lutetia. Recently, Vcrnazza ct al. [2009] 
and (partly) Nedelcu et al. [2007] showed that ECs are a 
good match for the visible/near-infrared spectra of Lutetia. 

The measured visual albedo for Lutetia has typically 
ranged over 15-22%, much higher than for the more 
common (CI/CM) carbonaceous chondrites (CC) and 
overlapping the lower range for S-types (in recent liter- 
ature, the early dedicated observations of Lutetia have 
been supplanted by reference to five rather inconsistent 
IRAS scans, which imply a still higher albedo and smaller 
effective diameter for Lutetia, well under 100 km, to which 
we assign less significance, especially because they are 
inconsistent with the mean size derived here). The radiom- 
etry by Mueller et al. [2006], reduced using two different 
thermal models, also yields albedos too high for most CC 
meteorites. A recent determination of visual albedo, using 
Hubble Space Telescope observations [Weaver et al. 2009] 
and the size/shape/pole results from the present paper 
and Carry et al. [2010b], indicate a value near 16%. This 
value, consistent with EC albedos (as well as metallic), 
is generally higher than most CCs, although some types 
of CCs, namely CO/CVs, have higher albedos, typically 
about 10%, with some COs getting as high as 15-17% 
[Clark et al. 2009]. 
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Radar observations of Lutetia (Magri et al. [1999, 
2007], confirnied by Shepard ct al. [2008]) showed that 
Lutetia has a moderate radar albedo (0.19-0.24), com- 
fortably in the mid-range of ECs, but considerably lower 
than metallic M-types and higher than most CCs. The 
uncertainty range in these values overlaps with some 
CO/CV carbonaceous chondrites compositions at the low 
extreme and with some metallic compositions at the high 
end. Hence, CO/CV carbonaceous chondrites cannot be 
ruled out based on albedo considerations alone. 

Indeed, numerous researchers in the last few years 
(Lazzarin et al. [2004, 2009, 2010]; Barucci et al [2005]; 
Birlanetal. [2006]; Barucci et al. [2008]; Perna et al. 
[2010], see summary by Belskaya ct al. [2010]) have argued 
that Lutetia shows certain spectral characteristics (e.g., 
in the thermal IR) that resemble CO and CV types and 
do not resemble a metallic meteorite. However, in these 
studies, the comparisons with EC meteorites was less thor- 
ough, partly because mid-infrared comparison data are not 
extensive. From rotationally resolved visible/near-infrared 
spectra of Lutetia, Ncdelcu et al. [2007] claimed a better 
match with CC in one hemisphere and with EC in the 
other, although this has yet to be confirmed. If Lutetia 
were highly heterogenous, that might explain some of the 
conflicting measurements. 

Vcrnazza ct al. [2010], however, have shown that 
mid-infrared emission of asteroids of similar composition 
can be very different due to differences in surface particle 
size. Miner alogical interpretations from this wavelength 
range are thus subject to caution and must be supported 
by VNIR reflectance spectra. Also, although some CO 
meteorites show albedos approaching that of Lutetia, the 
lack of a 1 micron olivine band in Lutetia's reflectance 
spectrum [see Fig. 3 of Barucci et al. 2005] argues against 
CO composition since most, but not all, COs have a 1 
micron band. Since the strength of this band generally 
shows a positive correlation with albedo, Lutetia's high 
albedo suggests that a strong 1 micron band should be 
present if its composition were CO. The lack of a drop-off 
in Lutetia's spectral reflectance below 0.55 micron and its 
high albedo make it inconsistent with CV meteorites [see 
GafFey 1976, for instance]. Finally, M and W-type asteroids 
(parts of the X class [DeMeo et al. 2009] if albedo is not 
known) have colors in the visible that are inconsistent with 
C-types. 

Colors, spectra, polarization, and albedos give us a 
picture of the relatively thin surface layers of an asteroid. 
Effects such as space weathering, repeated impacts that 
churn the regolith, recent impacts that may locally expose 
fresh material, particles sizes, or even differentiation 
processes may hinder our ability to ascertain the bulk 
composition of an object. Bulk density, on the other hand, 
gives us a picture of the entire asteroid body and ought 
to be a powerful constraint on bulk composition (subject 
to uncertainties about porosity and interior structure, 
mentioned below). Our new size estimates, when combined 
with recent mass determinations from other workers, 
now allow us to make estimates of the bulk density for 
Lutetia. Table ■') lists the volumes from three of the models 
addressed in this work, the triaxial ellipsoid model, the 
KOALA radius vector model, and our best estimate 
hybrid model. When coupled with two mass estimates, by 
Baer et al. [2008] or Ficnga et al. [2009] (which themselves 
differ by 25%), we find the given bulk densities. 



Table 5. Lutetia Mass, Volume and Density 



Method 


Vol 
(xlO^^cm^) 


Density (g cm '') 
P P 


Triax 
KOALA 
Hybrid 


6.46 ±0.95 
5.13 ±1.02 
5.94 ± 0.90 


3.98 ±0.69 3.19 ±1.04 
5.00 ±1.10 4.01 ±1.41 
4.32 ±0.77 3.46 ±1.13 



" with mass of 2.57 ± 0.24 x 10^^ g from Baer ct al. [2008] 
*> with mass of 2.06 ± 0.60 x 10^^ g from Fienga et al. [2009] 



Grain densities for stony meteorites range from 
^2.3 g.cm"'^ for CI/CM carbonaceous chondrites 
[Consolmagno et al. 2008] to -3.0-3.6 g.cm'^ for CO/CV 
[Flynn ct al. 1999] to ~3.6 g.cm'^ for EC [Mackc et al. 
2009]. Our best model yields densities of 4.3 or 3.5 g.cm^'^, 
which are among the higher densities yet tabulated for 
asteroids. The maximal extent of uncertainties on our 
preferred model range from about 2.3 to 5.1 g.cm"'^. 
Conservatively, if we were to consider meteorite grain 
densities, this range excludes iron-nickel, CI, and CM. 
Enstatite chondrites are favored, but CO/CV are also al- 
lowed. Perhaps more realistically, meteorite bulk densities 
should be considered instead. These are significantly lower 
- e.g., CO/CV appear to be, on average, about 0.6 g.cm~^ 
lower [Macke et al. 2009]. 

One must also be careful in trying to put too much em- 
phasis on comparison of asteroid densities with meteorite 
densities. Most asteroids are thought to have significant 
macroporosity [e.g., they may be rubble piles like (25143) 
Itokawa Fujiwara et al. 2006], so the asteroid density 
is likely to be substantially lower than the component 
material density [see Britt et al. 2006]. If so, then metallic 
meteorites are still ruled out, but only marginally at the 
upper end of the uncertainty range. CO/CV are similarly 
ruled out at the lower end of the uncertainty range, but 
not by much. The meteorite density values can only be 
guidelines. This provides some constraint on the possible 
bulk composition, but without reliable, smaller uncertain- 
ties in mass estimates, we must also rely on other observed 
quantities, such as albedo and spectra. A better mass 
estimate from Rosetta will reduce the density uncertainty 
considerably. 

In summary, our consideration of density and other 
evidence favors EC composition for Lutetia, and although 
CO/CV composition is not ruled out definitively, we con- 
sider it a lower probability. Among all the known meteorite 
classes, hydrated enstatite chondrites seem to fit the most 
number of measured parameters. These chondrites are 
represented among known meteorites only by the hydrated 
EC clasts in the unusual meteorite Kaidun. Finally, we 
emphasize that Lutetia may well be composed of material 
that is either rare or not yet represented in our meteorite 
collections. One example that might work is a low-albedo 
carbonaceous matrix material to suppress the olivine 
bands, embedded with abundant high-albedo clasts (such 
as an AUende-like composition, but with a much higher 
abundance of C AIs) . 



7. Summary 

We used adaptive optics images of (21) Lutetia from 
various large telescope facilities, and at various epochs, to 
make a triaxial ellipsoid shape model. In a companion pa- 
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per, we combine these AO images with lightcurves covering 
several decades to produce a radius vector model. There are 
advantages and disadvantages to these two methods. Here, 
we have combined the best properties of each to yield a 
hybrid shape model, approximated by an ellipsoid of size 
124 X 101 X 93 km (with uncertainties 5x4x13 km) that can 
be easily used to compute sizes, volumes, projected areas, 
and densities. When coupled with recent mass estimates, 
this hybrid model suggests a density of 3.5 ± 1.1 g cm^"^ 
or 4.3 ± 0.8 g cm~^. This is within the range expected for 
EC-like compositions, although the uncertainties formally 
permit other compositions. 

The Rosetta mission presents a unique opportunity for 
us to perform the ultimate calibration of our PBD and tri- 
axial ellipsoid approach to determine sizes and rotational 
poles. Furthermore, it will offer a chance to compare and 
contrast our triaxial ellipsoid model to the KOALA model 
for Lutetia. 
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